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ABSTRACT 

We use a wide- field (0.9 deg^) X-ray sample with optical and GALEX ultraviolet observations to 
measure the contribution of Active Galactic Nuclei (AGNs) to the ionizing flux as a function of redshift. 
Our analysis shows that the AGN contribution to the metagalactic ionizing background peaks at around 
z = 2. The measured values of the ionizing background from the AGNs are lower than previous estimates 
and confirm that ionization from AGNs is insufficient to maintain the observed ionization of the inter- 
galactic medium (IGM) at > 3. We show that only sources with broad lines in their optical spectra 
have detectable ionizing flux and that the ionizing flux seen in an AGN is not correlated with its X-ray 
color. We also use the GALEX observations of the GOODS-N region to place a 2a- upper limit of 0.008 
on the average ionization fraction /,y(700 A)//i,(1500 A) for 626 UV selected galaxies in the redshift 
range z = 0.9 — 1.4. We then use this limit to estimate an upper bound to the galaxy contribution in the 
redshift range z = — 5. If the z ~ 1.15 ionization fraction is appropriate for higher redshift galaxies, 
then contributions from the galaxy population are also too low to account for the IGM ionization at the 
highest redshifts (z > 4). 

Subject headings: cosmology: observations — cosmology: diffuse radiation — galaxies: active — 
galaxies: intergalactic medium 



general population. Bergvall et al. (2006) claimed a sig- 
nificant escape fraction in the local blue compact galaxy 
Haro 11. Shapley et al. (2006) detected ionizing fluxes 
from two z = 3 Lyman Break Galaxies (LBGs), though 
these measurements were at a much lower level than a pre- 
vious measurement made by Steidel et al. (2001). Iwata 
et al. (2008), using a narrow band filter technique, found 
ionizing radiation in seven z > 3 LBGs and 10 z > 3 Ly- 
man alpha emitters out of a sample of 198 z > 3 galaxies 
in the SSA22 field. However, their analysis of the level of 
contamination by foreground galaxies is quite crude, and 
it is also possible that at least some of the detections cor- 
respond to objects with misidentified redshifts. Paradoxi- 
cally, Iwata et al. found a highly significant null detection 
for the object SSA22a-D3, which is the brighter of the two 
objects for which Shapley et al. (2006) claimed an ionizing 
flux detection. 

Here we use UV images from the Galaxy Evolution Ex- 
plorer (GALEX) mission (Martin et al. 2005) to measure 
the actual ionizing fluxes from very large samples of both 
galaxies and AGNs. As we illustrate in Figure 1, the 
1528 A filter on GALEX samples the ionizing region of 
the spectrum at z = 1. The 2371 A filter samples it at 
z = 2. This avoids issues of the calibration of the spec- 
tra and also of the selection functions present in the sam- 
ples used to construct the composite spectra. It is also 
now recognized that complete samples of AGNs (at least 
those which are not Compton thick) are most easily ob- 
tained with X-ray observations. Barger et al. (2003) and 
Fontanot et al. (2007) have used X-ray samples to obtain 
direct limits on very high-redshift (z = 3.5 — 6.5) ionizing 
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Foundation. 
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1. INTRODUCTION 

One of the most important parameters in the cosmolog- 
ical modeling of galaxies and the intergalactic gas is the 
level of the ionizing radiation in the universe. This is often 
referred to as the metagalactic ionizing background. The 
metagalactic ionizing background is usually computed us- 
ing composite quasar spectra convolved with the quasar 
luminosity function (e.g., Haardt & Madau 1996; Madau 
et al. 1999; Haardt & Madau 2001; Meiksin 2005). Al- 
though the relative contributions of the galaxy popula- 
tion and the Active Galactic Nucleus (AGN) population 
are very poorly determined, based on the ionization stages 
seen in the intergalactic gas, we suspect that AGNs dom- 
inate the production below redshifts of about three. At 
higher redshifts recent estimates suggest that the AGNs 
are unable to account for the ionizing flux (e.g., Bolton et 
al. 2005; Meiksin 2005). In this paper we will confirm this 
result and show that the AGN contribution is even lower 
than previously calculated. 

It is usually assumed that the ionizing flux at these 
higher redshifts is predominantly produced by galaxies. 
However, there is very little evidence that enough ionizing 
photons escape from the galaxies for this to be true. Most 
measurements for local and z = 1 galaxies have given only 
upper limits on the escape fractions (Leitherer et al. 1995; 
Malkan et al. 2003; Siana et al. 2007), as have recent mea- 
surements at z = 3 (Giallongo et al. 2002; Inoue et al. 
2006; Fernandez-Soto et al. 2003). However, there is some 
counter evidence for significant escape fractions, though 
it is unclear how representative these galaxies are of the 
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Fig. 1. — The GALEX filter responses plotted for z = 1.15 (col- 
ored) and the intrinsic spectrum of a galaxy with a decaying star 
formation rate (black) having an exponentiation time of 5 X 10^ yr 
and an age equal to that of the universe a,t z = 1.15. The galaxy 
model is from Bruzual & Chariot (2003) and does not contain the ef- 
fects of the galaxy's neutral hydrogen opacity or internal extinction. 



In this paper we combine GALEX observations with X- 
ray samples to determine the ionizing flux from z ^ I 
AGNs and then use this result to estimate the AGN con- 
tribution to the metagalactic ionizing background in the 
redshift range z = 0—5. We also combine GALEX observa- 
tions with a large optical galaxy sample with spectroscopic 
redshifts to determine limits on the escape of ionizing pho- 
tons at z = 1.15. We then use this limit to constrain the 
galaxy contribution in the redshift range z = — 5. The 
structure of the paper is as follows. In §2 we describe 
our optical and X-ray samples and how we determined 
the ultraviolet fluxes. In §3 we examine the contribution 
of optically selected galaxies and X-ray selected AGNs to 
the ionizing flux. Only the X-ray selected AGNs are posi- 
tively detected. In §4 we determine how the ionizing flux 
relates to other properties of the X-ray selected AGNs, 
showing that the ionizing radiation is only seen in those 
with broad-lines in their optical spectra. In §5 we examine 
the evolution of the metagalactic ionizing flux. We show 
the ionizing fluxes as a function of redshift and compare 
our results to previous calculations. Wc summarize our 
results in §6. 

We assume Hm — 0.3, Q\ ~ 0.7, and Hq = 
70 km s~^ Mpc~^ throughout. All magnitudes are given 
in the AB magnitude system, where an AB magnitude is 
defined by ttiab = — 2.51og/^ — 48.60. Here fi, is the flux 
of the source in units of ergs cm"^ s~^ Hz"^. 

2. ULTRAVIOLET FLUXES FOR THE SAMPLES 

2.1. Optical Sample 

We take as our initial optical sample all F435W < 26 
galaxies in the 145 arcmin^ area of the ACS Great Ob- 
servatories Origins Deep Survey-North (GOODS-N; Gi- 
avalisco et al. 2004) field. The F435W, F606W, F775W, 
and F850LP magnitudes are taken from the ACS catalogs, 
the U magnitudes are from the [/-band images of Capak et 
al. (2004), and the spectroscopic redshifts are from Barger 
et al. (2008). Further details and catalogs may be found 
in Barger et al. (2008). 



The GALEX mission obtained a deep 150 ks exposure 
of the GOODS-N region in early 2004. We measured the 
near-ultraviolet (NUV; 2371 A central wavelength) and 
far-ultraviolet (FUV; 1528 A central wavelength) mag- 
nitudes at the positions of the optical sample using the 
GALEX images, which we obtained from the Multimis- 
sion Archive at STScI (MAST). Given the large point 
spread function of GALEX (4'.' 5 - 6" FWHM), we used 
an 8"' diameter aperture to measure the magnitudes us- 
ing the GALEX zeropoints of 20.08 for the NUV im- 
age and 18.82 for the FUV image from Morrissey et al. 
(2007). We measured the magnitudes for the brighter ob- 
jects (F435W = 20 - 23.5) in both 8" and 24" diameter 
apertures and used the median offset of —0.41 between 
these measurements to correct all of the 8" magnitudes to 
approximate total magnitudes. The offset is slightly larger 
than the —0.36 correction from 7.6" diameter magnitudes 
to total magnitudes determined by Morrissey et al. (2007). 
The magnitudes agree on average to within 0.05 mag with 
the SExtractor (Bertin & Arnouts 1996) magnitudes given 
in the GALEX NUV-f FUV merged catalog for the region 
in the GR4 data release. 

We measured the noise level in the images by randomly 
positioning apertures on blank regions of the sky and mea- 
suring the dispersion. (The procedure for selecting the 
blank regions is outlined in §3.1). We found la limits of 
26.8 in the NUV image and 27.4 in the FUV image. Of the 
6458 galaxies in the F435W sample, 1017 have FUV < 25.5 
and 2528 have NUV < 25 (with considerable overlap be- 
tween the two). 

However, because of the large point-spread function 
(PSF) of GALEX, contamination by neighbors is a seri- 
ous problem. This can be dealt with statistically, as we 
discuss in §3, but it is also of considerable interest to look 
at the properties of the individual galaxies. We therefore 
generated a subsample of isolated galaxies where we elimi- 
nated any source which was closer than 8" to a source with 
a brighter GALEX magnitude in the NUV band or where, 
based on a visual inspection, the position was clearly con- 
taminated by the wings of a nearby bright GALEX source. 
A substantial fraction of the galaxies in the optical sample 
are eliminated by this isolation requirement, reducing the 
effective area to 48 arcmin^ based on the fractional number 
of sources rejected by the isolation criteria. The isolated 
galaxy sample contains 966 sources with NUV < 25 and 
422 sources with FUV < 25.5. The spectroscopic identifi- 
cations of these sources are 96% complete to NUV = 24.5 
and 98% complete to FUV = 25. 

2.2. X-ray Samples 

We construct our X-ray samples from the 2 — 8 keV 
sources in the Chandra Large Area Synoptic X-ray Survey 
(CLASXS) of Yang et al. (2004), in the Chandra Lock- 
man Area North Survey (CLANS) of Trouille et al. (2008), 
and within 8' of the pointing center of the 2 Ms Chan- 
dra Deep Field-North (CDF-N) survey of Alexander et al. 
(2003). We take as our full X-ray sample the sources with 
/2-8 kcv > 3.5 X 10"-^^ ergs cm~^ s~^ (note that only 
the CDF-N contributes below 7 x 10~^^ ergs cm~^ s~^). 
We take as our bright X-ray sample the sources with 
/2-8 kcV > 7 X 10^^^ ergs cm~^ s~"^. Current tables of 
redshifts for the full CLASXS, CLANS, and CDF-N X-ray 
samples are given in Trouille et al. (2008). We define X- 
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ray AGNs within the sample on energetic grounds as any 
source more luminous than Lx ~ 10"*^ ergs s"^ (Zezas et 
al. 1998; Moran et al. 1999). Here Lx is the luminosity 
calculated in the rest-frame 2 — 8 keV band. 

The full X-ray sample contains 662 X-ray sources. Al- 
most all of the sources have been spectroscopically ob- 
served and 442 have spectroscopic redshifts. The spectro- 
scopically identified sources should contain a nearly com- 
plete sample of all optical broad-line AGNs in the region 
(Barger et al. 2005; Richards et al. 2005). We find 171 
broad-line AGNs, of which 98 have X-ray quasar luminosi- 
ties {Lx > 10^** ergs s^^). The bright X-ray sample with 
/2-8 kov > 7 X 10~^^ ergs cm~^ s~^ contains 485 sources, 
of which 390 have spectroscopic redshifts. Here we have 
162 broad-line AGNs, of which 97 have X-ray quasar lu- 
minosities. 

The CLANS and CLASXS fields are covered by a num- 
ber of GALEX pointings, all of which are roughly 30 ks in 
depth, somewhat shorter than the 150 ks CDF-N GALEX 
exposure. The la limits are slightly variable but typically 
about 26.1 in the NUV images and 26.5 in the FUV images. 
Only a small fraction (about 4%) of the X-ray sources in 
the CLANS or CLASXS fields are not covered by one of 
the GALEX pointings and all of the CDF-N sources are 
covered. The effective area for the X-ray sources lying 
within the GALEX pointings is 0.9 deg^ at /2-8 kcV above 
7 X 10~^^ ergs cm^^ s^^ and 0.055 deg^ below. 

We measured the GALEX magnitudes for the sources 
in the X-ray samples which lay within the GALEX point- 
ings in the same way that we measured the GALEX mag- 
nitudes for the GOODS-N optical sample. Because the 
X-ray sources are generally brighter in the FUV than the 
optical sources, the problem of contamination by galaxies 
is smaller, and we did not eliminate sources using a separa- 
tion criterion. (The X-ray sources themselves are so sparse 
that contamination by another X-ray source may be ne- 
glected.) However, we did eliminate four sources where vi- 
sual inspection showed that there was contamination from 
a nearby brighter GALEX source at the same wavelength. 
None of these four sources appear to have any significant 
FUV flux centered on the object position. 

2.3. Redshift- Magnitude Relations 

In Figure 2 we show redshift versus FUV magnitude for 
(a) the isolated optical sample with FUV < 25 and (b) 
the bright X-ray sample with FUV < 25. In both cases 
nearly all of the sources are spectroscopically identified, as 
can be seen from the incompleteness histograms shown at 
the bottom of each panel. We have redshifts for 323 of the 
329 sources with FUV < 25 in the optical sample. Five of 
the remaining sources could not be identified, and the final 
source has not been observed. All of the 164 sources with 
FUV < 25 in the bright X-ray sample have been observed, 
and only 5 are not identified. In both panels we denote 
X-ray AGNs by red solid squares, and we enclose X-ray 
quasars in red large open squares. 

At z = 0.6 (dashed horizontal lines in Figure 2) the 
FUV filter straddles the Lyman break and substantially 
samples the ionizing radiation (see Fig. 1). As we move to 
higher redshifts nearly all galaxies vanish from the sam- 
ple. We can see from Figure 2a that the only source at 
z > 0.7 in the optical sample with a significant ionizing 
flux is a broad-line quasar. There is one additional galaxy 



at z > 0.7, but the FUV flux in this object is likely from 
a nearby z = 0.316 galaxy, which lies 2.7" away. (This 
source was not eliminated from the sample because the 
companion is fainter in the NUV band.) 
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Fig. 2. — Redshift vs. FUV magnitude for (a) the isolated op- 
tical sample (GOODS-N) and (b) the bright X-ray sample (X-ray 
sources in the CLANS, CLASXS, and CDF-N fields with /2_8 keV > 
7 X 10~^^ ergs cm~^ Sources with Lx > 10'*^ ergs s~^ 

(X-ray AGNs) are denoted by red solid squares, and those with 
Lx > 10** ergs (X-ray quasars) are shown enclosed in red 
large open squares. The fraction of spectroscopically unidentified 
sources per 0.5 mag bin is shown in histogram form at the base of 
each panel. The z = 0.6 redshift at which the FUV filter straddles 
the Lyman break is shown by a dashed horizontal line in each panel. 



We show this in a different way in Figure 3, where 
we plot the redshift distributions of all the spectrosopi- 
cally identified sources in the F435W < 24.5 GOODS- 
N sample (black histogram) and of only those detected 
at FUV < 25.5 (blue shading). (The GOODS-N sample 
is nearly spectroscopically complete to the F435W=24.5 
limit.) We also plot the redshift distributions of all the X- 
ray AGNs present in the optically selected F435W < 24.5 
GOODS-N sample (green histogram) and of only those de- 
tected at FUV < 25.5 (red shading). Roughly 76% of the 
galaxies in this sample below z = 0.6 have FUV < 25.5. 
Of the 378 sources in this sample above z = 0.7, only 
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four galaxies and two AGNs, including the one galaxy and 
the one broad-line quasar which we previously noted, have 
FUV < 25.5. 
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Fig. 3. — The black histogram shows the spectroscopic redshift 
distribution of the F435W < 24.5 sample in the GOODS-N field. 
The binning interval is 0.05 in redshift. Galaxies which are detected 
with FUV magnitudes less than 25.5 are shown with the blue shaded 
histogram. The green histogram shows the spectroscopic redshift 
distribution of the X-ray AGNs in this sample. X-ray AGNs which 
are detected with FUV magnitudes less than 25.5 are shown with 
the red shaded histogram. 
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Fig. 4. — Spectral energy distributions of the seven broad-line 
quasars with Lx > 10** ergs s~^ in the GOODS-N region hav- 
ing redshifts in the range z = 0.8 — 3. Each quasar is shown with 
a unique symbol. The GALEX points are shown in purple. The 
24 ^ra fluxes are from Treister et al. (2006) and are shown in red. 
The purple solid line shows the Lyman break, and the purple dashed 
line shows the position of Lya. The black solid line shows an fu 
vs. v~'^-^ power law. Only one source has a significant ionizing flux 
(open squares), and it appears to have almost no absorption. 

The FUV-bright broad-line quasar is only one of a num- 
ber of broad-line quasars in the GOODS-N region, and it 
is immediately clear that these quasars are extremely var- 
ied in the number of ionizing photons that we see from 
them. In Figure 4 we show the spectral energy distri- 
butions (SEDs) of the seven broad-line quasars with red- 
shifts between z = 0.8 and z = 3 in the GOODS-N region. 
Nearly all of the quasars have significant breaks across the 



ionization edge, while the one quasar where we see a signif- 
icant ionizing flux (open squares) shows a smooth exten- 
sion from its longer wavelength SED and almost no opacity 
to the ionizing photons. A fraction of the breaks will be 
produced by intervening Lyman Limit Systems (LLS). We 
quantify this in §4. 

The bright X-ray sample in Figure 2b, with its much 
larger area than the GOODS-N optical sample, contains a 
significant number of sources with detected ionizing fluxes 
that we can use to analyze the distribution in properties 
of the AGN ionizers and their contribution to the ionizing 
background. 

3. CONTRIBUTIONS OF OPTICALLY SELECTED GALAXIES 
AND X-RAY SELECTED AGNS 
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Fig. 5. — (a) Distribution of FUV fluxes in ^Jy for the isolated 
optical sample in the GOODS-N with z = 0.9 — 1.4 (weak or no X- 
ray emission - black histogram; AGNs - red histogram). The smooth 
blue curve shows the distribution of fluxes measured at random po- 
sitions in the FUV image. All but the one broad-line quasar, which 
has a 19 /ijy flux, are shown in the plot, (b) Same distribution for 
our bright X-ray sample with z = 0.9 — 1.4 (red histogram). Note 
the much wider x-axis. 

In Figure 5a we show the FUV flux distribution of 
the isolated optical sample in the GOODS-N with z = 
0.9 — 1.4. In Figure 5b we show the FUV flux distribu- 
tion of the bright X-ray sample with z = 0.9 — 1.4. The 
lower redshift cut-off of z = 0.9 places the FUV filter well 
below the Lyman break. We have chosen a conservative 
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upper limit on the redshift of z = 1.4 to place the filter 
well above the drop-off in the intrinsic galaxy spectrum at 
about 500 A (see Fig. 1) and to minimize the eflfects of 
intervening intergalactic absorption. 

3.1. Optical Galaxy Sample 

In Figure 5a we show the FUV flux distributions of the 
sources with weak or no X-ray emission (black histogram) 
and of the X-ray AGNs (red histogram) in the isolated 
optical sample in the GOODS-N with z = 0.9 - 1.4. We 
overplot the distribution of fluxes measured at random iso- 
lated positions in the GOODS-N FUV image (blue curve), 
which matches extremely well to the black histogram. The 
one detected galaxy in Figure 2a is seen at 0.3 /iJy, but 
the one detected broad-line quasar is far off the sc;ale of 
the plot. Excluding only the broad-line quasar while keep- 
ing the other X-ray AGNs, these isolated galaxies formally 
contain a total positive signal of 3.0 ± 0.8 /iJy. However, 
when we assess this result we must be concerned about 
contamination by foreground galaxies and systematic ef- 
fects of the background subtraction. Indeed, the average 
FUV flux per galaxy is only 0.011 /xJy, which is consider- 
ably below the la flux of 0.039 /iJy in the image. 

These effects are best tested by Monte-Carlo simula- 
tions using random position samples. This also allows us 
to avoid any subjectivity that might be present in the iso- 
lation criteria, since we can analyze analytically selected 
samples and treat the effects of blending using the simu- 
lations. The present problem is also simplifled since the 
contaminating effects primarily arise from the z ^ 0.6 — 0.7 
galaxies which, as we have seen in Figure 3, are by far the 
most common class of source with substantial FUV emis- 
sion. Since we do not expect these to be correlated with 
the 2; ~ 1 galaxies of interest, we do not need to be con- 
cerned with the effects that might be introduced by corre- 
lations, and we can use spatially uniform random samples 
to determine our backgrounds. 

The simplest procedure is to create an analytically se- 
lected isolated galaxy sample. To do this we masked 
16" radius regions around all z = — 0.7 galaxies and 
stars with FUV < 23 and 8" radius regions around all 
z = — 0.7 galaxies and stars with FUV < 25.5 and ex- 
cluded all z = 0.9 — 1.4 galaxies which lay in the masked 
regions from our sample. In total, we masked regions 
around 469 galaxies and 14 stars. We also masked 16" ra- 
dius regions around the four X-ray AGNs in the field with 
FUV < 25.5. These isolation criteria removed 187 galaxies 
from our initial sample of 680 galaxies with 2; = 0.9 — 1.4 
in the GOODS-N region and reduced the effective area to 
104 arcmin^. The average FUV flux per galaxy in the ana- 
lytic isolated sample is 0.012 /iJy, which is similar to what 
we found for the smaller isolated sample we constructed 
in §2.1. 

We next created 100 random realizations, each contain- 
ing the same number of positions as the target sample. We 
constrained the random positions to lie away from the low- 
redshift galaxies, stars, and AGNs in the same way that 
we constrained the target sample. However, since the sig- 
nal associated with the z = 0.9 — 1.4 galaxies is small, we 
did not remove or avoid the actual sample when construct- 
ing the random samples, but simply generated random x 
and y positions within the unmasked area. From the 100 
random realizations we found that a random position sat- 



isfying the isolation criteria has an average FUV flux of 
0.018 /iJy and a dispersion of 0.003 /iJy. Subtracting this 
from the average FUV flux per galaxy of 0.012 fj. Jy mea- 
sured from the analytic isolated target sample gives an 
FUV flux of -0.006 ± 0.003 ^Jy per galaxy. 

The average NUV flux per galaxy in the same sample 
is 0.44 /iJy, while the average J7-band flux per galaxy is 
0.96 fxJy. (The corresponding fluxes for a random position 
at these wavelengths are negligible compared to these val- 
ues.) The rest-frame wavelengths of the two bands at the 
midpoint of the z = 0.9 — 1.4 redshift interval are 1100 A 
and 1720 A, respectively. The NUV band straddles the 
Lyman continuum break in the higher redshift portion of 
the redshift interval which will reduce the observed flux 
in this bandpass and it may also be reduced by the Ly- 
man alpha forest (see Fig. 1). We therefore chose to use 
the f/-band flux as a measure of the rest-frame flux at 
1500 A, assuming a flat fi, SED. We then formed the ra- 
tio of the rest-frame 700 A flux to the rest-frame 1500 A 
flux, ignoring the small differential if -correction as a func- 
tion of redshift. We call this ratio the ionization fraction 
to distinguish it from the escape fraction of ionizing pho- 
tons from the galaxy. Conversion of the ionization fraction 
to the escape fraction requires knowledge of the intrinsic 
galaxy spectrum and reddening (e.g., Haardt et al. 1999; 
Siana et al. 2007). However, for the present purposes of 
obtaining the ionizing emission of the galaxy population, 
it is the ionization fraction which is of interest and which 
we can use to convert the rest-frame UV emission to the 
ionizing emission. The ionization fraction we obtained us- 
ing the background-corrected analytic isolated sample is 
-0.006 ±0.003. 

An alternative way to proceed is to subtract the fore- 
ground galaxies and stars to form a cleaned image and then 
to measure the signal associated with the target sample on 
the cleaned image. This has the virtue of allowing us to 
measure a larger fraction of the target sample and to avoid 
using the somewhat arbitrarily sized masking regions se- 
lected in the previous procedure. To form a cleaned image 
we fitted and subtracted all of the z < 0.7 galaxies with 
FUV < 25.5 in the GOODS-N from the GALEX image, 
as well as all of the stars and the four X-ray AGNs with 
FUV < 25.5. As can be seen from Figures 2a and 3, this 
should remove nearly all of the signal from directly de- 
tected objects and therefore minimize the effects of blend- 
ing and optimize the background subtraction. 

Our procedure for this cleaning was to normalize the 
GALEX PSF, which we determined by stacking the FUV 
detected stars in the image, to the measured flux of 
each galaxy and then to subtract this from the image. 
Most galaxies are very compact compared to the broad 
GALEX PSF, but to remove the 36 brighter objects with 
FUV < 22.5, which may be more extended, we masked 
out a 12" region around these galaxies. In our subsequent 
analysis we exclude objects which fall within these masked 
regions and correct the area appropriately. 

For each subset sample from the full F435W < 26 sam- 
ple to be analyzed, we measured the fluxes in the cleaned 
image at the actual positions of the sample and then cre- 
ated 100 random samples with an equal number of po- 
sitions as the sample to determine the background and 
dispersion. We then subtracted this background from the 
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measured signal of the sample to obtain the true signal. 
Again, since the signal associated with the z — 0.9 — 1.4 
galaxies is small, we did not remove or avoid the actual 
sample when constructing the random samples. 

For the 626 spectroscopically identified z = 0.9 — 1.4 
galaxies with F435W < 26 that are not X-ray AGNs and 
are not in the small masked areas, we find an average 
background-corrected FUV flux per galaxy of —0.005 ± 
0.002 /iJy. The average NUV flux per galaxy is 0.51 /iJy, 
and the average [/-band flux per galaxy is 0.97 /j,Jy. We 
note again that the NUV and U signal associated with 
a random postion is negligible. The inferred ionization 
fraction is —0.006 ± 0.003. This procedure gives a con- 
sistent result with the analytic isolation simulations de- 
scribed above. Both give negative fluxes but are consistent 
with a zero flux at just over the 2tT level. 

We can also directly compare the distribution of fluxes 
in the observed target sample with that from the total 
of all 100 realizations of the random sample to test if 
they are different from one another. In Figure 6 we com- 
pare the cumulative distribution of the FUV fluxes for the 
z = 0.9 — 1.4 F435W < 26 sample measured from the 
cleaned image (black curve) with that measured from the 
random sample (red curve). A Kolmogorov-Smirnov test 
give a 1% probability that the target sample has the same 
distribution as the random sample. 
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Fig. 6. — Cumulative fraction of sources vs. FUV flux measured 
in the cleaned image. The black curve shows the F435W < 26 
sources in the z = 0.9 — 1.4 redshift interval, and the red curve 
shows the random sample. 

For the nearly spectroscopically complete F435W < 
24.5 sample, there are 363 z — 0.9 — 1.4 galaxies that 
are not X-ray AGNs and are not in the small masked ar- 
eas. Here we find an average background-corrected FUV 
flux per galaxy of -0.005 ± 0.003 /.tJy. The average NUV 
flux per galaxy is 0.65 /iJy, and the average C/-band flux 
per galaxy is 1.34 /iJy. The inferred ionization fraction is 
-0.005 ±0.003. 

We should probably not take the negative signal too se- 
riously, given its statistical significance in all of the mea- 
surements described above. However, the average nega- 
tive flux may seem worrying at first sight and suggestive 
of some type of systematic error in the procedure. We 



tested this by measuring the signal for the 229 galaxies 
with z > 1.6 with the analytic isolation procedure. Here 
we find an average background-corrected FUV flux per 
galaxy of 0.001 ± 0.004 /iJy. Thus, the negative signal 
seems to be associated only with the specific target sam- 
ple. 

However, it is worth noting that the galaxies could have 
an effective negative signal if their shadowing effect on the 
ionizing background exceeded their direct emission. The 
galaxies will block any background light at these wave- 
lengths over the region where their neutral hydrogen col- 
umn density makes them a LLS. The shadowing effect may 
extend over a larger area than any FUV emission, making 
the larger apertures needed for the GALEX observations 
more sensitive to it. Such shadow galaxies could be used 
to make a measure of the ionizing flux and of the extent 
of the neutral hydrogen in the individual galaxy. How- 
ever, higher spatial resolution observations would be more 
appropriate for such an analysis than the GALEX data. 

In all of our estimates above, systematics could raise the 
absolute value of the error, though they will not change the 
ratio of the signal to the error. In particular, we should 
allow for uncertainties in the extrapolation to form the 
1500 A flux, possible systematic errors in the determina- 
tion of the FUV fluxes (Morrissey et al. 2007) , and uncer- 
tainties in the aperture corrections to total magnitudes. 
Specifically, we include the possibility of a systematic er- 
ror of 0.2 in the relative FUV and C/-band magnitude de- 
terminations, which would raise the error in the ionization 
fraction to 0.004. Allowing for this level of systematic er- 
ror and also ignoring the negative signal, we take the 2a 
limit of 0.008 as our upper limit on the ionization fraction 
in the z ^ 1 galaxy population in what follows. 

Given the relatively small size of the GOODS-N field, 
we might also be concerned that we are missing the rare 
FUV-bright galaxies. Since we see no FUV-bright galaxies 
in the field, we adopt a la upper limit of 1.8 such galax- 
ies (Gehrels 1986). Adopting a typical value for the FUV 
magnitude of 23 — comparable to the brightest z ~ 0.6 
galaxies (Fig. 2) — one such bright galaxy would corre- 
spond to an average contribution of 0.006 /iJy per galaxy 
from the 680 measured galaxies. This is already smaller 
than the adopted limit. However, using an FUV magni- 
tude of 23 is in fact extreme, since even young star-forming 
galaxies will have a strong break across the Lyman edge, 
so the true contribution will be smaller. It does remain 
conceptually possible that there may be unique environ- 
ments where UV ionizing galaxies preferentially form and 
that we simply have not sampled these with the GOODS- 
N area. This can only be tested with larger area samples. 

In Figure 7a we show the ratio of the rest-frame 700 A 
flux to the rest-frame 1500 A flux (the ionization frac- 
tion) as a function of the absolute rest-frame 2000 A 
magnitude computed from observed-frame F435W (black 
squares). Here we use the cleaned FUV images and do 
the background subtraction as described above. The red 
line shows the 2a upper limit of 0.008 on the average ion- 
ization fraction for the full F435W < 26 sample. In Fig- 
ure 7b we show the ionization fraction versus observed- 
frame F435W - F850LP color (roughly rest-frame 2000 A 
- 4200 A). There is no significant detection as a function 
of either luminosity or color. 
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Fig. 7. — (a) Average ionizing fraction (ratio of the rest- frame 
700 A flux to the rest-frame 1500 A flux) vs. absolute rest-frame 
2000 A magnitude computed from observed-frame F435W. The 
black squares show the results for the spectroscopically identified 
z = 0.9 — 1.4 galaxies in the F435W < 26 sample that are not X- 
ray AGNs. The error bars are Ic statistical errors. The red solid 
line shows the 2it upper limit of 0.008 on the average ionization 
fraction for the full sample, including possible systematic errors de- 
scribed in the text, (b) Average ionizing fraction vs. observed frame 
F435W-F850LP color (roughly rest-frame 2000 A - 4200 A) for the 
z = 0.9 — 1.4 galaxies. 



3.2. Bright X-ray Sample 

The one FUV-bright broad-line quasar in the GOODS- 
N gives a flux of 19 /iJy, which dominates the signal in 
the field. We show this in Figure 8, where we compare 
the summed image formed using our cleaned image of the 
626 spectroscopically identified z — 0.9 — 1.4 galaxies in 
the F435W < 26 sample that are not X-ray AGNs (left 
panel) with the summed image of the 29 GOODS-N X-ray 
AGNs with z = 0.9 — 1.4 (right panel). The clear detection 
in the FUV sum of the X-ray AGNs is totally dominated 
by the FUV-bright broad-line quasar. The remaining 28 
X-ray AGNs in this redshift interval are consistent with 
having no ionizing flux. However, given the small number 



of contributing broad-line AGNs, we need the much wider 
field of the bright X-ray sample to accurately compute the 
AGN contribution to the ionizing flux. 

In Figure 5b we show the FUV flux distribution for our 
bright X-ray sample with z = 0.9 — 1.4 (red histogram). 
32 sources in this redshift interval have FUV magnitudes 
brighter than the 5a limit of the image, and the total mea- 
sured flux is 112 fiJy for the 0.9 deg^ area. This result is 
not affected by our exclusion of the four sources which 
are contaminated by neighbor galaxies (see §2.2), since, 
even without any cleaning, these sources only have a total 
flux of 2 fi,Jy. We analyzed the background level and vari- 
ance by again creating random realizations, each contain- 
ing the same number of positions as the target sample, and 
measuring at the positions in the raw images. We found 
the background to be small relative to the signal. The 
background-corrected signal is 100 ± 5 /iJy, where the er- 
ror is the variance in the background determination. This 
measured signal would correspond to 4.5 /iJy in the small 
GOODS-N region. Thus, the one FUV-bright broad-line 
quasar in the GOODS-N region is an anomaly. 

Even with the larger area, there are only a small num- 
ber of AGNs which are seen to have ionizing radiation. 
The three brightest AGNs produce 60% of the flux, and 
the seven brightest AGNs produce almost 80%. Given the 
small number of contributing sources and the peculiar dis- 
tribution of the fluxes, we have used a jackknife analysis 
to make a more accurate estimate of the uncertainty in the 
AGN contribution. This gives a 68% confidence range of 
±33 ^Jy on the signal. It should be noted that large un- 
certainty in the jackknife analysis does not mean that the 
ionizing flux contribution from the AGN is not highly sig- 
nificant. This is properly inferred from the results in the 
last paragraph. The large uncertainty from the jackknife 
method is a consequence of the fact that only a small frac- 
tion of the X-ray sources which are detected in the FUV 
dominate the production of ionizing radiation. 

We conclude that only the X-ray sources can be de- 
tected at ionizing frequencies and that an extremely small 
number of these dominate the ionizing flux production. 
The emission per unit comoving volume of the X-ray 

-1) is uX^ = (3.4 ± 1.5) X 
: 1.15, where we have used 



quasars {Lx > 10 ergs s 



IQ-i^ ergs s~ Mpc~'^ at z 
the jackknife method to estimate the uncertainty. Here Aj/ 
is the luminosity density per unit frequency, which is the 
sum of the luminosities per unit frequency of all the ob- 
jects divided by the cosmological volume that they occupy. 
We now turn to a more detailed analysis of the contribu- 
tion from the bright X-ray sample to try to understand 
the FUV-bright AGNs in more detail. 

4. IONIZING PROPERTIES OF X-RAY SELECTED AGNS 

We can try to refine the results of §3 by considering how 
the ionizing flux relates to other properties of the AGNs. 
In Figure 9a we show Lyi/(700 A), where 700 A is the 
rest- frame wavelength, for the z — 0.9 — 1.4 sources in the 
full X-ray sample versus rest-frame 2 — 8 keV luminosity, 
Lx- We divide the sources into those with broad lines (red 
squares) and those without (black diamonds). As might 
be expected, only the broad-line AGNs are ionizers. One 
of the non-broad-line AGNs is significantly detected, but 
an inspection of its spectrum shows that it appears to be 
contaminated by a lower redshift galaxy. Thus, we exclude 
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Fig. 8. — Summed image of the 626 spectroscopically identified z = 0.9 — 1.4 galaxies in the F435W < 26 sample that are not X-ray AGNs 
and do not lie in the small number of masked areas in the cleaned image (left panel) and the 29 GOODS-N X-ray AGNs with z = 0.9 — 1.4 
(right panel). The noise is higher in the galaxy image because of the larger number of objects. The scaling is the same in both images. In 
each case the image size is 80" on a side, and a positive signal would appear at the center of the image. The strong detection in the right 
panel is dominated by the one FUV-bright broad-line quasar. There is no detected signal in the left panel. 



it from further consideration. 

In Figure 9a wc use large open squares to show the mean 
values of L,yi'{700 A) versus Lx for the z = 0.9 — 1.4 broad- 
line AGNs. For all of the z = 0.9 - 1.4 broad-hne AGNs 
together we find a mean ratio of X,^z^(700 A)/Lx = 3.95. 
However, there is a very wide spread in the individual 
ratios. The dominant sources have ratios that arc al- 
most an order of magnitude higher than the mean ratio 
of 3.95. The FUV and optical are more tightly corre- 
lated (sec Fig. 9b). so a considerable part of the observed 
spread in Figure 9a comes from the dispersion in the rela- 
tion of the optical to the X-ray fluxes of the sources (see 
Fig. 9c). A linear relation with a normalization of 3.95 
(blue dashed line) does not describe the mean values of 
-L^i^(700 A) with Lx very well. However, a non-linear 
relation does. The best power law fit to the mean val- 
ues gives L^i^{700 k)/Lx = 3.53 x {Lx/IQ'^'^ ergs s-i)"-^^ 
with a la uncertainty of 0.05 in the power law index (black 
solid line). This non-linear dependence follows from the 
well-known non-linear dependences of the optical and ul- 
traviolet luminosities of broad-line AGNs on X-ray lumi- 
nosity (e.g., Vignali et al. 2003a, 2003b; Strateva et al. 
2005; Steffen et al. 2006; see Fig 9c). 

In Figure 9b we show the dependence of Lj,v{7Q0 A) 
on i^i/(2300 A) for the z = Q.Q - 1.4 sources in the full 
X-ray sample, which is tighter. The rest-frame 2300 A 
luminosity is computed from the (/-band magnitudes ob- 
tained on the CFHT Mcgaprime camera by Trouille et 
al. (2008). This filter is centered at an observed- frame 
wavelength of 4900 A. Once again we ignore small differ- 
ential X-corrcctions assuming is flat. Wc sec it is the 
broad-line AGNs that make up the UV-luminous popu- 
lation. The broad-line AGNs divide about equally into 
those which are well described by a ~ j^-i-SS pQ-^^gr 
law between the two wavelengths (green dotted line) and 
those which have little or no ionizing flux. Wc use large 
open squares to show the mean values of Lvv{7Q0 A) ver- 
sus L^v{2mQ A) for the .2 = 0.9 - 1.4 broad-line AGNs. 
For all of the z = 0.9 — 1.4 broad-line AGNs together we 
find a mean ratio of vL^{100 k)/vL^{2m) A) = 0.47. A 
linear relation with this normalization (blue dashed line) 
provides a reasonable fit to the mean values. In fact, the 
best power law fit to the mean values (black line), which 
has a power law index of 0.90 ±0.18, is consistent with the 



linear relation. 

We can compare our mean ratio with that of Trammell 
et al. (2007), who analyzed the GALEX fluxes for a large 
sample of Sloan Digital Sky Survey (SDSS) quasars. Their 
quasars are about an order of magnitude higher in lu- 
minosity than the present sample, but their mean ratio 
of i/L^ (700 A)/z/L^(2300 A) is similar. Trammell et al. 
(2007) only included sources with detected FUV magni- 
tudes, for which they find a ratio of about 0.5. Allowing 
for the roughly one-third of sources which they do not de- 
tect, this value could be as low as 0.33. However, given 
that these two numbers closely bracket our present results, 
it is clear that there is not a large luminosity dependence 
in the ratio. 

Finally, in Figure 9c we show the well-known non-linear 
dependence of Ljxi^(2300 A) on Lx- The best fit power 
law gives vL^{2?,m)/Lx = 7.46 x {Lx/IQ'^^ ergs ^-^f-^^ 
with a 1(7 uncertainty of 0.21 in the index (black line). 
This is almost identical to the relation found by Vignali 
et al. (2003a; sec also Richards et al. 2005 for details on 
converting the Vignali et al. relation into the form used 
here) (red dotted line). 

One can next ask if there is any relation between the 
X-ray properties of the broad-line AGNs and the ioniz- 
ing photon escape. In Figure 10 we plot the ratio of 
L^i>{7QQ k)/Lx for > 3 x 10''^ ergs s"^ sources with 
redshifts z = 0.9 — 1.4 versus X-ray softness, as measured 
by the ratio of the 0.5 - 2 keV flux to the 2 - 8 keV 
flux. We divide the sources into those with broad lines 
(red squares) and those without (black diamonds). We 
can again see that all of the non-broad-line AGNs are UV 
faint. This immediately tells us that we are not misidenti- 
fying sources as non-broad-linc AGNs when they are really 
broad-line AGNs, as one might have expected to happen 
if the broad lines were not visible spectroscopically due to 
dilution by the host galaxy (Moran et al. 2002; Severgnini 
et al. 2003; Cardamone et al. 2007). Barger et al. (2005) 
make a similar argument based on the presence of nuclei 
in the broad-line AGNs. 
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Fig. 9. — The z = 0.9 — 1.4 full X-ray sample (broad-line AGNs - red solid squares; non-broad-line AGNs - black diamonds; moans for the 
broad-line AGNs - red open squares), (a) Li^u at rest-frame 700 A, Li,u{70Q A), vs. rest-frame 2 — 8 keV luminosity, Lx- (b) L^i/(700 A) 
vs. iji/i-'{2300 A). The green dotted line shows the relation for a source with ~ (j.) L^i^(2300 A) vs. Lx- The red dotted line 

shows the relation converted from Vignali et al. (2003a). In all three panels we show linear relations (blue dashed lines) and best power law 
fits (diagonal black lines) (see text). In (a) and (b) the 3a errors in L^j/(700 A) for the CLANS and CLASXS fields are shown by the dashed 
horizontal lines. In (a) and (c) the Lx above which the field size corresponds to the wide area is shown by vertical dotted lines. Below this 
the sources are only in the much smaller area of the CDF-N. Sources with low or negative y-axis values are plotted at a nominal value of 
1.2 X 10''^ ergs s'^ . The only non-broad-linc AGN with significant y-axis values in (a) and (b) appears to be contaminated by a lower redshift 
galaxy (see text). 
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Fig. 10. — Ratio of L,^u at rest-frame 700 A, Lyi'(700 A), to 
rest-frame 2 — 8 keV luminosity, Lx, vs. X-ray softness, as mea- 
sured by the ratio of the 0.5 — 2 keV flux to the 2 — 8 keV flux, 
/o.5-2 kev//2-8 keV for the z = 0.9 - 1.4 redshift range. Only 
sources with 2 — 8 keV luminosities above 3 X 10'*'^ ergs s~^ are 
shown. Broad-line AGNs are denoted by red squares, and non- 
broad-line AGNs by black diamonds. The flux ratio for a power law 
spectrum with photon index F = 1.8 is shown by the vertical dashed 
line. Sources with 700 A fluxes less than 0.2 fijj are shown with 
downward pointing arrows. 

Indeed, we can pick out a large fraction of the broad- 
line AGNs simply on the basis of their having a signifi- 
cant FUV flux. Specifically, none of the 45 non-broad-line 
AGNs with Lx > 10*^ ergs s'^ and redshifts z = 0.9 - 1.4 
are detected above the FUV flux limit of 0.26 /xJy (Sct), 
while 70 ± 14% of the broad-hne AGNs are. Thus, the 
GALEX observations can be used to pick out most of the 
broad-line AGNs in the X-ray sample. More interestingly, 
the selection does not seem to depend on the X-ray soft- 
ness, though the number of X-ray hard broad-line AGNs 
in the sample is small. As can be seen in Figure 10, we 
have six broad-line AGNs with 0.5-2 keV to 2-8 keV flux 
ratios less than 0.3, and four of these have significant FUV 
detections. There seems to be nothing unusual about the 
optical spectra of the broad-line AGNs which are detected. 
Apparently, broad-line AGNs can have a high ratio of the 
ionizing flux to the X-ray flux, even if the source is quite 
X-ray hard. This at once says that there is not a one-to- 
one correspondence between the X-ray spectral index and 
the neutral hydrogen opacity to the source, since any sub- 
stantial neutral hydrogen opacity would have the effect of 
wiping out the ionizing flux, as indeed is invariably seen 
in the non-broad-line AGNs. A detailed discussion of the 
comparison between the X-ray and the optical and ultavi- 
olet properties will be given in L. Trouille et al. (2008, in 
preparation). 

Regardless of the interpretation of the results, we con- 
clude that the production of the ionizing flux is best traced 
by optical spectroscopic broad-line AGN selection, rather 
than by X-ray or combined X-ray and optical selection, 
since the X-ray color of the broad-line AGN is not cor- 
related with the presence of the ionizing flux, and non- 
broad-line AGNs are simply not detected. This consider- 
ably simplifles the task of computing the evolution of the 
metagalactic flux with redshift, since we believe that the 



spectroscopic broad-line AGNs are relatively completely 
identified in the X-ray samples. 

The remaining question is why some broad-line AGNs 
are essentially unobscured in the FUV (we shall refer to 
these as clear-channel broad-line AGNs), while others with 
similar luminosities are substantially absorbed in the FUV 
(we shall refer to these as extinguished broad- line AGNs). 
A substantial part of the explanation is obscuration by 
intervening LLSs along the line of sight which, even at 
z = 1, is a significant effect. At z = 1 the density of LLS 
is approximately 0.9 per unit redshift, which would result 
in about 40% of the sources being obscured (e.g. Storrie- 
Lombardi et al. 1994). This is very similar to the fraction 
of sources with weak ionizing flux in the present sample 
and in Trammell et al. (2007). 

However, a detailed analysis shows the situation is 
slightly more complicated. Storrie-Lombardi et al. (1994) 
flnd that Mg II absorption line systems with rest-frame 
equivalent widths above 0.3 A are a very close proxy for 
LLS. Thus, we can approach the problem empirically, since 
we can easily detect such systems in the spectra of our 
broad-line AGNs. 

In Figures 11a and lib we show the spectral index a 
(where ~ ly") computed between 700 A and 2300 A 
versus L^iy{2300 A). In Fi gure 11a we use the observed 
FUV flux to compute the rest-frame 700 A flux, and in 
Figure lib we use the observed NUV flux. We do not in- 
clude the additional /iT-correction as a function of redshift 
relative to the mean redshift, since this would require as- 
sumptions about the SED. However, this effect should be 
small. In Figure 11a we show the broad-line AGNs with 
z = 0.9-1.6andL^j/(2300 A) > 5 x 10^^^ ergs s'^. In Fig- 
ure lib we show the broad-line AGNs with z = 1.9 — 2.5 
and L^iy{2300 A) > 2 x 10^^ ergs s'^. Note that the up- 
per redshift limit in the latter is set by the requirement 
that Mg II at the quasar redshift be on the spectrum. The 
spectral indices of sources with FUV fluxes below the 2a 
threshold have been calculated using the 2a threshold and 
are shown with downward pointing arrows. At lower lu- 
minosities the uncertainties in the FUV fluxes become too 
large to determine accurately the spectral indices. We use 
red squares to show sources without Mg II absorbers and 
we use black squares to show sources with Mg II absorbers 
that are associated with the AGNs. We denote sources 
with intervening Mg II absorbers along the line of sight 
by black solid diamonds when the redshift of the system is 
such that it would substantially quench the FUV light and 
by black open diamonds when the redshift of the system is 
such that it only partially covers the wavelengths observed 
by the FUV filter. 

As expected, the presence of an intervening Mg II ab- 
sorber at a redshift which would substantially quench the 
FUV light (the black diamonds) guarantees that a source 
is extinguished. (We take sources with spectral indices 
less than —2 to be extinguished.) However, two of the 
clearly extinguished sources in the lower redshift interval 
(Fig. 11a) do not have Mg II absorbers. In one case the 
wavelength coverage of the spectrum does not extend to 
short enough wavelengths to ensure that there is not a 
LLS along the line of sight that could be significantly ex- 
tinguishing the flux in the FUV fllter. In the other case 
there is no Mg II absorption, and if the extinction is caused 
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by a LLS, then it must have weak Mg II absorption. The 
possible very small fraction of LLS without Mg II absorp- 
tion is discussed in Storrie-Lombardie et al. (1994) and 
Steidel & Sargent (1992). 
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Fig. 11. — (a) Spectral index computed between 700 A and 
2300 A for 2 = 0.9 - 1.6 broad-line AGNs with rest-frame 2300 A 
luminosities above 5 X 10** ergs s~^ vs. L,^u at rest-frame 2300 A, 
(2300 A). Sources with 700 A fiuxes below the 2a threshold are 
shown at the 2a threshold with downward pointing arrows. Sources 
without Mg II absorbers in the line of sight are shown with red small 
squares. Sources with Mg II absorbers associated with the AGNs 
are shown with black larger squares. Sources with intervening Mg II 
absorbers along the line of sight are shown with black solid diamonds 
when the rcdshift of the system is such that it would substantially 
quench the FUV light and with black open diamonds when the red- 
shift of the system is such that it would only partially cover the 
wavelengths observed by the FUV filter. The black line shows the 
median spectral index for all of the sources. The blue dashed line 
shows the median spectral index for the sources after removing those 
with spectral indices less than —2 that do not have Mg II absorbers 
associated with the AGNs. (b) Same plot for z = 1.9—2.5 broad-line 
AGNs with rest-frame 2300 A luminosities above 2 X 10*^ ergs . 

Interestingly, for four of the sources with Mg II ab- 
sorbers in Figure 11a, the Mg II system is a narrow- line 
system associated with the AGN (within 10,000 km s~^ of 
the AGN) and is not intervening (the four black squares). 



In three of the cases the velocity separation is less than 
2000 km s^^ from the associated AGN. For one other 
source in Figure lib (the black square) the quasar shows 
broad absorption features. All of these cases show some 
level of absorption but in two cases we still see an ionizing 
flux. Either the conversion of Mgll to neutral hydrogen 
in these systems is different from that of the intergalac- 
tic sytems or they only partially cover the ionizing region. 
The actual explanation is not important for the present 
purposes since there is a critical distinction between as- 
sociated and intervening systems. In the associated cases 
the ionizing radiation is absorbed before it emerges into 
the general intergalactic medium. In contrast, intervening 
systems are part of the radiative transfer in the intergalac- 
tic medium. Thus, in computing the ionizing source term 
we should not correct for absorption by the associated sys- 
tems while we should remove the effects of absorption in 
the intervening systems. 

The median spectral index of all of the sources in the 
lower rcdshift interval is shown by the black line in Fig- 
ure 11a and has a value of —1.8 (-2.7,-1.4), where the 
bracketed indices give the 68% confidence range computed 
using the median sign method. If we remove the three sys- 
tems with substantial intervening absorption (three black 
solid diamonds) and the two clearly extinguished sources 
without Mg II absorption (two red squares)^ assuming they 
are also caused by an intervening LLS, then the median 
spectral index becomes —1.45 (—1.9,-1.2) (blue dashed 
line) and the mean spectral index is —1.35. The num- 
ber of sources in the higher rcdshift interval (Fig. lib) is 
small, but the median spectral index after removing the 
two systems with substantial intervening absorption (two 
black solid diamonds) is consistent with that in the lower 
rcdshift interval. Thus, we assume the mean intrinsic spec- 
tral index is rcdshift invariant and adopt the mean value 
of —1.35 at z = 1.15 in computing the evolution of the 
metagalactic ionizing flux in §5. Our 700 A to 2300 A flux 
ratio is 0.20, and we scale this to a 912 A to 2300 A flux 
ratio of 0.26, assuming an intrinsic spectrum with an in- 
dex of — 1 below the break, which is similar to that of the 
broad-line AGNs above the break. This is a factor of 1.5 
lower than the value adopted in Haardt & Madau (1996). 
Using our value in their analysis would bring their results 
closer to the present values. 

5. EVOLUTION OF THE METAGALACTIC IONIZING FLUX 

The discussion of the previous sections suggests that 
the best procedure to measure the metagalactic flux is to 
take the broad-line AGN sample from the X-ray data and 
to measure its light at NUV wavelengths above the Ly- 
man break. We can then use the average spectral break, 
corrected for intergalactic absorption, to derive the ioniz- 
ing flux cmissivity. Our broad-line AGN sample chosen 
from the full X-ray sample probes to much fainter UV 
magnitudes than the large optical sample and therefore 
avoids major extrapolations in determining the UV vol- 
ume emissivity at high redshifts. Using the measured UV 
fluxes rather than normalizing to the X-ray fluxes avoids 
the scatter and non-linearity in the relation of the UV to 
the X-ray. Finally, measuring the flux at the longer wave- 
lengths and then correcting to the ionizing flux avoids the 
complex problem of correcting for the intervening inter- 
galactic absorption, though at the expense of assuming 
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a redshift invariant average spectrum for the broad-line 
AGNs. 

In Figure 12 we show the evolution of the comoving 
volume emissivity determined in this way from the broad- 
line AGNs in the CLANS, CLASXS and CDF-N fields (red 
squares). We first computed the rest-frame 2300 A emis- 
sivity for our broad-line AGNs in each redshift interval. 
To determine the rest -frame 2300 A flux, we interpolated 
(or extrapolated at the the highest redshifts) using the 
NUV magnitudes and the g (observed frame 4900 A) and 
i (observed-frame 7900 A) magnitudes from Trouille et al. 
(2008). We then scaled the 2300 A emissivity to the ion- 
izing volume emissivity using the results from §4. The 
errors were again computed with the jackknife method. In 
order to test the completeness, we computed the results as 
a function of the limiting X-ray flux. Below z = 2.5 the 
results are dominated by sources with X-ray fluxes above 
7 X 10^"'^^ ergs cm^^ s^^, which are sampled by the larger 
area. Above this redshift the results are dominated by 
sources in the smaller CDF-N area and cosmic variance 
could be a more significant issue. However, the results 
substantially converge at fluxes well above the limiting flux 
in the CDF-N sample. We do not detect any broad-line 
AGNs above z — A, implying a large fall-off in the ionizing 
fluxes above z — -i. The evolution of our comoving volume 
emissivity can be parameterized by the equation 

A„„ = (1.67-|-15.91z-3.660z2)1023 ergs Hz^^Mpc^^ , 

(1) 

which is shown as the solid red curve in Figure 12. 
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Fig. 12. — Redshift evolution of the comoving volume emis- 
sivity of the ionizing flux just below the Lyman continuum edge 
from broad-line AGNs. The red squares show the values calculated 
directly from the broad-line AGNs in the full X-ray sample (see 
text). The error bars are Icr calculated using the jackknife method. 
There are no broad-line AGNs above z = 4 in our sample, and we 
show the upper limit based on assuming a one sigma upper limit of 
1.8 sources (Gehrels 1986) with the characteristic break luminosity 
shown in Figure 13b. The blue curve was calculated from the max- 
imum likelihood fit to the Richards et al. (2005) optically selected 
broad-line AGN sample. We show this only over its valid redshift 
range. The black solid curve shows the source function computed by 
Meiksin (2005), which was estimated from the SDSS survey data for 
broad-line AGNs combined with lower luminosity broad-line AGN 
constraints from an optical survey at z = 3. We have corrected this 
to the presently derived spectral shapes. 



In the above calculation we have worked directly with 
the 2300 A emissivities determined from the broad-line 
AGN sample in the CLANS, CLASXS, and CDF-N fields. 
A similar result can be derived from broad-line AGN X-ray 
luminosity functions determined as a function of redshift 
(e.g., Yencho et al. 2008 and references therein), but this 
would require the use of the more uncertain conversion 
from X-ray emissivity to ionizing emissivity. 

We can also determine the results from the luminosity 
functions of optically selected samples of broad-line AGNs, 
where these are deep enough. Here we have computed 
the comoving ionizing volume emissivity of the broad-line 
AGNs from the maximum likelihood fit to the Richards et 
al. (2005) optically selected sample (blue curve). As would 
be expected from the fact that Richards et al. (2005) found 
their luminosity functions to be in almost perfect agree- 
ment with the X-ray luminosity functions, the results de- 
rived from the optical and X-ray samples are in extremely 
good agreement over z = — 2, where the optical sam- 
ples are sufficiently deep to make this comparison. It is 
much harder to use the existing wide-field optically se- 
lected samples of broad-line AGNs to make such a com- 
parison at high redshifts, since such samples are not deep 
enough to probe the dominant contributors. Extrapola- 
tions outside the range for which a fit has been derived 
can produce wildly inaccurate results as discussed by for 
example Richards et al. (2006). 

Meiksin (2005) used the SDSS observations of broad- 
line AGNs together with the lower luminosity broad-line 
AGN constraints from the optical survey of Steidel et al. 
(2002) and Hunt et al. (2004) at z = 3 to estimate the 
higher redshift values from the optical samples. We show 
his results with the black solid curve. (We have modified 
the Meiksin 2005 result which assumed a slightly different 
value of the break across the Lyman edge to the presently 
derived spectral shape.) Given the extrapolation uncer- 
tainties in this type of estimate, the agreement is quite 
good. As is well known, the results of Haardt & Madau 
(1996) and Madau et al. (1999) are too high because of 
the assumed quasar luminosity functions. However, the 
present results are considerably lower than even the modi- 
fied Haardt-Madau values that are often used (e.g., Haardt 
& Madau 2001; Bohon et al. 2005). 

As we have discussed, the production of the ionizing 
flux is dominated by a small number of sources. We can 
investigate this directly using the present X-ray selected 
sample. Typically 3—10 sources produce most of the flux 
in our line of sight in the redshift intervals of Figure 12, 
even with a 1 deg^ area. (Different lines of sight would 
see different AGNs as being broad-line and having signif- 
icant ionizing flux because of the geometric effects of the 
intrinsic absorption in the AGNs. However, the bulk of the 
luminous X-ray sources are broad- line AGNs [e.g., Steffen 
et al. 2003; Barger et al. 2005], and so we expect this 
effect to be small.) We quantify the number of contrib- 
utors in Figure 13a, where we show the number density 
of sources producing 67% of the observed ionizing flux. 
The number density drops from a value of 10~^ Mpc^'^ 
at low redshifts, where lower luminosity sources dominate, 
to about 10^^ Mpc""^ a.t z — 3.5. In Figure 13b we show 
the luminosity range that produces from 25 — 75% of the 
ionizing flux. This emphasizes that a rather narrow range 
of luminosities straddling the break luminosity in the lu- 
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minosity function produces most of the flux. The typical 
luminosity of the dominant sources rises quite rapidly at 
low redshifts, reflecting the rapid nearly pure luminosity 
evolution of the optical and X-ray luminosity functions in 
this redshift range (e.g., Barger et al. 2005; Richards et 
al. 2005). It then reaches a fairly invariant characteris- 
tic luminosity of about 7 x 10^^ ergs s^^ Hz^^ above this 
redshift. This corresponds to an absolute rest-frame mag- 
nitude at 2300 A of M2300 = -24.5. 
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Fig. 13. — (a) Redshift evolution of the comoving volume density 
of the sources producing 67% of the ionizing volume emissivity. The 
1(7 error bars have been computed using the median sign method. 
The number of dominant sources in each redshift interval ranges 
from 3 to 10. (b) The luminosity range of the sources dominating 
the ionizing volume emissivity. The squares show the luminosity 
above which 50% of the emissivity is produced, and the vertical 
lines show the 25% to 75% range. 

As pointed out by Madau et al. (1999), the measured 
absorption distance for the ionizing photons at z ~ 3 is 
quite small so that the ionization rate is only determined 
by local sources. They estimate an absorption distance 
corresponding to Az — 2.8(1 -t- at these redshifts, 
which is Az = 0.18 at z = 3. Combining this with the 
comoving number density in Figure 13 gives a contribut- 
ing number of quasars of about 3000 at z = 3 and 200 
at z = 4 in any given region. Beyond this the numbers 
drop rapidly, and at z = 4.5 there would be less than 20 
dominant quasars. We may compare this with Meiksin 
& White (2003), who self-consistently estimated the ab- 
sorption distance by normalizing the modeled quasars to 



match the required ionization in the Lya forest and then 
computing the absorption distance, rather than obtaining 
it from observations. The number of contributing quasars 
is higher by factors of several, as would be expected since 
the present AGN emissivity is too low to account for the 
ionization at these redshifts. As the number of contribut- 
ing quasars drops, fluctuations become important, and the 
contribution of the AGNs to producing transmission in the 
AGNs is further reduced, increasing the deficiency of these 
sources for ionizing the IGM (Meiksin & White 2003). 
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Fig. 14. — Redshift evolution of the comoving volume emis- 
sivity of the ionizing flux just below the Lyman continuum break 
calculated from the broad-line AGNs in the full X-ray sample (red 
squares) compared with the upper limits for galaxies (blue and green 
squares). The blue solid square shows our directly measured 2cr 
limit ai z = 1.15. The green open squares show the 2a limits ob- 
tained by converting the 1500 A comoving emissivities from Tresse 
et al. (2007) to the ionizing emissivities using our 2(t upper limit 
on the average ionization fraction at z = 1.15 (see §3.1) and then 
multiplying by 1.26 to correct from 700 A to 912 A. The black dia- 
mond shows the measured value from Shapley et al. (2006)'s LBGs. 
The dotted line shows the required emissivity to reproduce the Lyo 
forest structure in the IGM estimated by Meiksin (2005). 

In Figure 14 we compare the comoving volume emissiv- 
ity of the ionizing light just below the Lyman continuum 
edge calculated from the broad-line AGNs in the full X- 
ray sample ( red squares ) with the upper limits for galaxies 
(blue and green squares). At z = 1.15 we show the directly 
measured 2a upper limit from the present work as the blue 
square with the downward pointing arrow. To compute the 
limits at higher redshifts, we converted the 1500 A comov- 
ing emissivities from Tresse et al. (2007) to the ionizing 
emissivities using our 2fj upper limit on the average ion- 
ization fraction at z = 1.15 (see §3.1) and then multiplied 
by 1.26 to correct from 700 A to 912 A. For the highest 
redshift points (z ^ 3) we used the unweighted estimated 
values from Tresse et al. (2007) 's Table 4. We show these 
limits on the comoving ionizing emissivity derived from 
the Tresse et al. (2007) data as green open squares with 
downward pointing arrows in Figure 14. The Tresse et 
al. (2007) 1500 A emissivities are broadly consistent with 
previous measurements (Steidel et al. 1999; Gabasch et 
al. 2004; Arnouts et al. 2005; Sawicki & Thompson 2006) 
and with the present work. The black dotted curve in 
Figure 14 shows the the required emissivity to reproduce 
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the Lyo; forest structure in the intergalactic medium at 
z > 3 estimated by Meiksin (2005). We have used the 
value allowing for reprocessing of the ionizing radiation in 
the IGM (the dashed curve of Meiksin's Figure 1). 

If the high-redshift ionization fraction is as low as our 
measured hmit at z = 1.15, then the upper hmit on the 
galaxy contribution is sUghtly low compared with the re- 
quired ionization. This result is marginal in the z = 3 — A 
range but stronger above z = 4. However, given the rapid 
change in the morphologies of the galaxy population be- 
yond z = 1, we may question whether the z — 1.15 ion- 
ization fraction of the present work is applicable to the 
higher redshift population. If the ionization fraction rises 
to that measured by Shapley et al. (2006) from their small 
sample of LBGs, then the galaxy contribution at high red- 
shifts, as shown by the black diamond at z = 3, would 
easily match (or indeed exceed) the required ionization. 
(We note that if Iwata et al. [2008] are correct in their null 
detection of the brighter of the two Shapley et al. [2006] 
LBG detections, then the Shapley et al. point would fall 
by a factor of roughly 3.) However, the Shapley et al. 
(2006) sample may be biased towards the LBGs that are 
most likely to have ionizing flux and thus may really be an 
upper limit even on the LBG contribution. It is therefore 
critical to understand whether the current measurements 
of the z ~ 3 LBGs are representative of the population 
as a whole at these redshifts. The narrow-band ionizing 
flux surveys of z 3 LBGs currently Ixung worked on by 
several groups should resolve this question. 

6. SUMMARY 

We have used X-ray, optical, and GALEX observations 
to measure the contribution of AGNs to the ionizing flux 
as a function of redshift. Our analysis of a large popu- 
lation of X-ray sources confirms that the AGN contribu- 
tion to the ionizing flux peaks at around z = 2 and drops 
rapidly at higher redshifts. It is insufficient to account for 



the observationally inferred ionizing flux at high redshifts 
(z > 3). We have also obtained a strong upper limit on 
the contribution of galaxies to the ionizing flux at z = 1.15 
using GALEX observations of a very large sample of rest- 
frame UV selected galaxies in the GOODS-N region. Our 
2(7 upper limit on the ionization fraction for this popula- 
tion, /;.(700 A)//^(1500 A) = 0.008, yields an upper hmit 
on the comoving ionizing emissivity from the galaxies at 
z = 1.15, which is at most comparable to that for the 
AGNs at z = 1.15. We briefly discuss the possibility that 
the galaxies may in fact be net absorbers at these redshifts. 

If galaxies are to contribute signiflcantly to the ioniz- 
ing radiation, then the ionization fraction must increase 
at higher redshifts. Shapley et al. (2006) 's measured value 
from their small sample of LBGs at z = 3 is about six 
times higher than the value from our 2a upper limit on 
the ionization faction at ^; = 1.15. However, the Shapley 
et al. (2006) value would be reduced by about a factor of 
three if the object which is not conflrmed by Iwata et al. 
(2008) is removed. If these LBGs are representative; of the 
high-redshift population as a whole, then this would be 
adequate to produce the high-redshift ionizing flux. How- 
ever, it would require an increase in the escape of the ion- 
izing photons from the galaxies as a function of increasing 
redshift. 
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